BNL-HET-03/14 



Lattice Matrix Elements and CP Violation in B and K 
Physics: Status and Outlook* 

Amarjit Soni 
High Energy Theory Group, Physics Department 
Brookhaven National Laboratory, Upton, NY 11973 
(email: soni@bnl.gov) 

ABSTRACT 

Status of lattice calculations of hadron matrix elements along with 
CP violation in B and in K systems is reviewed. Lattice has pro- 
vided useful input which, in conjunction with experimenatl data, 
leads to the conclusion that CP-odd phase in the CKM matrix plays 
the dominant role in the observed asymmetry in B —> ipKs- It is 
now quite likely that any beyond the SM, CP-odd, phase will cause 
only small deviations in B-physics. Search for the effects of the new 
phase(s) will consequently require very large data samples as well as 
very precise theoretical predictions. Clean determination of all the 
angles of the unitarity triangle therefore becomes essential. In this 
regard B — > KD° processes play a unique role. Regarding K-decays, 
remarkable progress made by theory with regard to maintenance of 
chiral symmetry on the lattice is briefly discussed. First application 
already provide quantitaive information on Bk and the AI =1/2 
rule. The enhancement in ReAo appears to arise solely from tree 
operators, esp. Q2', penguin contribution to RcAq appears to be 
very small. However, improved calculations are necessary for e'/e as 
there the contributions of QCD penguins and electroweak penguins 
largely seem to cancel. There are good reasons, though, to believe 
that these cancellations will not survive improvements that are now 
underway. Importance of determining the unitarity triangle purely 
from K-decays is also emphasized. 

I Introduction 

With important input from the lattice along with the classic results of indirect 
CP violation in Kl — > tttt, the asymmetric B-factories with measurements of CP 
asymmetry in _B ^ 4^Ks are providing valuable support to the CKM paradigm 
of CP violation[Tj. It is now clear that the CP-odd phase in the CKM matrix is 
the dominant source of CP violation in B ipK^- However, as is well known 
essentially compelling theoretical arguements suggest that new CP-odd phase(s) 

"Invited talk at 9th International Symposium on Particles, Strings and Cosmology (PAS- 
COS 03), Mumbai (Bombay) India, 3-8 Jan 2003 



should exist due to physics beyond the SM (BSM). At the same time there is no 
good reason to think that their effects in B-physics would be particularly large. 
Indeed, SM teaches a valuable lesson in this regard: even though the CKM phase 
causes a huge asymmetry (i.e. 0(1)) in i? ^ V'-^^s, its effects in CP violation 
in K-decays is miniscule « 10~^. Clearly this realization should motivate us to 
prepare for small deviations from the predictions of the SM in B-Physics even 
if the new CP-odd phase is large. For this reason, not only we need very large 
data samples of B's giving impetus to Super-B factories along with BTeV and 
LHCB, we also need exteremely precise tests of the SM. Residual theory errors 
are a serious cause of concern as they can easily thwart experimental efforts for 
search of BSM CP-odd phase(s)|g. 

With this perspective in mind, a brief discussion of the lattice method and 
results for the hadronc matrix elements that are important for weak interaction 
phenomenology are given in Sections I to IV. Therein I also discuss some of the 
anticipated experimental input that could help attain greater precision in con- 
straining CKM parameters. Section V emphasizes concern about residual errors 
in theory and the importance of clean determinations, i.e. without theoretical 
assumptions, of all the angles of the unitarity triangle. In this regard the special 
role of i? ^ KD'^ processes is also emphasized there. 

Progress made in the past few years with regard to maintenance of exact 
chiral symmetry on the lattice is outlined in Section VI along with application 
of this development to Bk- 

Section VII gives a brief report on the results from the 1st application of 
domain wall fermion method, which exhibit excellent chiral behavior, to K ^ 
TTTT, A/ —1/2 rule and e'/e. These 1st applications give good insight to the 
A/ —1/2 rule; in particular, contribution of penguin operators to ReAo in our 
lattice calculations appears to be extremely small and most of the enhancement 
seen in RcAq is originating from the tree operator, Q2. Unlike in the case of 
the A/ = 1/2 rule, the approximation currently used though appear too crude 
to give reliable information on e'/e. This difficulty arises as contributions of 
QCD penguins and electroweak penguins substantially cancel. However, there 
are very good reasons to suspect that this cancellation is not "natural" and is 
unlikely to survive as calculations are improved. 

For the purpose of stringent tests of the CKM model of CP violation a 
separate determination of the unitarity triangle purely from K-decays, to be 
compared to that obtained from B-physics, is highly desirable and this is finally 
emphasized in Section VIII. 

II Lattice Methodology: a very brief recapitu- 
lation 

Recall, Green's functions are calculated by numerical evaluation of the Feynman 
path integral. 
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- / WiQ)Mdet M{U)exp[~Sg{U)] 
^ ^ / DU det MiU)exp[-Sg{U)] ^ ' 

As it stands, dependence of the quark matrix M in this expression on the Unk 
variables renders its evaluation extremely difficult. To fecilitate the numerical 
calculation, one often uses the quench approximation (QA)and sets detM = 1. 
Physically, this approximation corresponds to neglect of the qq vacuum polar- 
ization loops in the propagation of the gauge field. The hint that this may be a 
reasonable approximation originally came from deep inelastic scattering exper- 
iments wherein the effect of qq pairs in the "sea" is accurate to about 15% [3]. 
There are, though, very good reasons that tell us that the accuracy of the QA 
in lattice computations is process dependent. 

In the past several years more and more "unqucnchcd" simulations, s i.e. 
those not using the QA so that dynamical qq pairs are included, have been 
underway. These studies show that QA seems to be valid to about 5 — 10% 
accuracy in non-singlet hadron spectrum 

On the other hand, dynamical quarks seem to increase the B-meson pseu- 
doscalar decay constant quite appreciably, (at least when rup is used to set the 
scale) [5] 

Ib'^^/Ib'^'^ = 1.23 ± 0.04 ± 0.06 (2) 

In addition to the QA there are several other sources of systematic errors 
in a typical lattice gauge calculation. Chief among these are finite (box) size 
and finite lattice spacing (a) errors. Also most lattice simulations are done with 
rather large values of masses of light (u,d,s) quarks and rather low values for 
the b quark mass, compared to their physical values. Painstaking and elaborate 
efforts become necessary to accurately extract from the data information rele- 
vant to the physical case. This may, for example, require extrapolation of the 
data (at a fixed gauge coupling, or lattice spacing) as a function of quark mass 
to the chiral limit and also extrapolation of the data as a function of the lattice 
spacing to the continuum limit (i.e. lattice spacing goes to zero). Furthermore, 
simulations for a fixed gauge coupling at two or more volumes are often needed 
for extrapolation to the infinite volume limit. 



II. 1 Some Examples of Brute-Force 

Relevant to this talk there are three works which serve to illustrate compu- 
tational brute force used in bringing them to fruition; these are Bk, fs and 
K- > 2tt. 

1. Bk- a major accomplishment of the lattice gauge effort, and in particular 
of the JLQCD group is their result0C|, Bk = 0.860 ± 0.058 in the QA. 
During the first 6-7 years ('84-'91 ), several exploratory attempts were 
made[S]. Methodology was in place around '91 [H] which was followed by 
several years of intensive computations leading to the final result obtained 
around '98. In the past few years this important result has been the focus 
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Figure 1: Bk vs. a, the lattice spacing. Only the 2 data points (circles) from 
DWF use non-perturbative renormalization of lattice operators, all others use 
one-loop perturbation theory. Of the circled DWF 2 points, the coarser spacing 
(right) corresponds to data from NERSC whereas the finer one (left) belongs to 
RBC collaboration [T^ 



of further checks and confirmation using other fermion discretizations, 
i.e. Wilson ^U] as well as the newer discretizations; domain wall fermions 
|11[I12|[T^ and the overlap fermions |14[I15| . The results of these methods 
are in rough agreement with the JLQCD result; however, with domain wall 
quarks (DWQ) method the central value of Bk tends to be 10-15% below 
the JLQCD result which may amount to a discrepancy of around 1-2 
a. (see Fig. ^ More precise calculations with these newer discretizations 
including a study with dynamical domain walls jl6| is now underway. 



2. /b- Another example is provided by fs, the B- meson pseudoscalar decay 
constant wherein the "heavy" b-quark mass « 4.5 GeV represents an ad- 
ditional technical problem. After the initial 5-6 years of exploratory works 
the computational strategy became quite well known around '92 '17'. In- 
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deed the result in the QA has been quite stable and withstood checks with 
the use of different techniques EI ■ In the past few years there has 
been some weak indication from experiment that heavy-light decay con- 
stant are somewhat smaller in the QA compared to experiment '2U' (i.e. 
full QCD). Indeed after years of persistent study the MILC collaboration 
has now finished the calculation of /b in quenched as well as in full QCD 
(i.e. with three light flavors of dynamical quarks) |^ and find, the 
ratio given in eq.(2). 

giving fs — 207 ± 35 MeV[22]. Clearly many independent checks and 
confirmations will take place in the next few years. 

3. -ftT — > TTTT and e'/e. Our third example is the calculation of the matrix 
elements of iiT — > 27r and e'/e using domain wall quarks (DWQ) by the 
CP-PACS [2S|and RBC collaborations [I3|. Unlike Bk and fs this is a 
first attempt to address if ^ 27r by both the collaborations in which not 
only QA but also a few other key approximations are made [see below]. 
Nevertheless, given the complexity of the problem it must be considered 
an inportant accomplishment which even at this early stage is providing 
very useful information on the long standing issue of the A/ = 1/2 rule. 
However, its repercussions on e'/e require careful study of systematic errors 
and improved calculations, which could take another few years. 

Ill Lattice Matrix Elements and CKM Constraints 

In the Wolfenstein representation, the CKM matrix can be parameterized in 
terms of the four parameters, A, A, p and fj 24 . Of these A = sin 6'c = 0.221 ± 
0.002, is the best known, A is known with modest accuracy, A — 0.847 ± 041 
and p and fj are poorly known. An important objective where lattice can help 
is in determination of fj and p accurately, fj is intimately related to the CKM 
phase i5i3 indeed SM cannot accommodate any CP violation if 77 = 0. 

The basic strategy is very simple. Assuming the SM is correct, and using the 
necessary theoretical input one translates experimental results on to an allowed 
domain on the fj - p plane. If a (new) experimental result requires value of p 
and/or fj that are inconsistent with those extracted from existing experiments 
then that could mean a failure of the SM. 

III.l Theoretical background and brief comments 

For the past several years, the following four experimental measurements have 
been used for an extraction of fj and p: 

1. The indirect CP violation parameter, e = (2.274 ± .017) x 10"^ 

2. The Bd - Bd mass difference, Ams^ = 0.487 ± 0.014^5"! 
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3. The Bs — Bg mass difference, for which at the moment only a lower bound 
exists, AniB, > IS.Ops"^ at 95% CL. 

This important bound is provided by experiments at LEP and SLD |26| . 
It is widely anticipated that an actual measurement of Ato^^ (rather 
than just a bound) will be accomplished at the Tevatron in the next few 
years. This will be very important for CKM determinations as the ratio 
Ams^/Ams^ can give Vtd/Vts if the SU(3) breaking ratio of hadronic 
mixing matrix element could be determined [27| . 

4. Ruc = b^ ulv/b -> civ = 0.085 ± 0.017. 
RecallEHl 

|e| = MKCKA^\mi^^S{xc) + mS{xt)A^\\l - p) + r^^S{x,, xt)} (3) 

Here Xq = m^/M^, where q = u,c,t i.e. the virtual quarks in the box graph 
for K'^ — K'^ oscillations, and S{xq) are the so called Inami-Lin functions |29|. 
Also, 

6V2n^AmK 

Mk - \mUK°\[si,{1 - l,)dYm (5) 

is essentially the hadronic matrix element. Once the Mk is known fi,p can be 
constrained through the use of eq.3. This matrix element is often parametrized 
in terms of Bk which should equal 1 if vacuum saturation approximation (VSA) 
holds. Since Jk (and rriK) is known quite precisely from experiment, evaluation 
of the matrix element is completely equivalent here to that of Bx- 
Similarly, we note that for Bd — Bd oscillation 

XBd = MBd CBd[{l -p)^+n ^]llQCDS{xt) A^\^/TBd (6) 

where XBd = p"'"'' and CBd — 5— f — §— Again once the hadronic matrix ele- 

ment, Mb^, is known eq. © can be used to constrain p, fj. 

This matrix element is a 3-point function, which is directly calculable on 
the lattice. More often than not, though, in analogy with the kaon case, M^d is 
parametrized in terms of a "B-parameter" defined as 

<Bd\\h,,[l~-1^)dY\Bd> 

^Bd - 8^2 f2 

3,"''BdJ Bd 

Then the physical quantity XBd requires both Bsd and fBd since the latter is not 
yet known from experiment. Besides since MBd seems to scale roughly as /^^ 
one needs to know fBd rather accurately. Also, in practice in most calculations 
of fBd one tries to fit the light quark mass dependence through some linear 
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function; such a fit, though is unlikely to give precisely the dependence on light 
quark mass for fg^. 

Once Bs — Bg oscillation are experimentally detected and AmB, becomes 
known, then the ratio 

Ama, ^ \Vta\'' {Ba\[b^,{l^^,)df\Ba) 

can be used to determine | Vtd\ if the ratio of hadronic matrix elements could 
be determined from the lattice. 

Since this ratio of matrix elements is completely dependent on SU(3) break- 
ing effects (s ^ d) it is expected to be close to unity. The objective of lattice 
calculations should be a precise evaluations of this SU(3) breaking and this 
necessitates an accurate treatment of light quarks. 

Again, introducing "B-parameters" for Bd and B^ mesons we can rewrite 

— = —^^m-pr+fi'] (9) 

where ^ is the SU(3) breaking ratio, 



^ = T^J^ (10) 

jBd V -DSd 

Finally, the semi-leptonic branching ratio b ulv/b — > civ 

is another important way of constraining p, 77 as it is a function of Vub/Vcb, 

.A(p-2 + f)/(l-^) (11) 



I Vcb\ ' 2 

To deduce Vub/Vcb, from the experimental measurement of the branching 
ratios requires corresponding form factors for exclusive reactions wherein lattice 
methods can be useful ^30, . In the interest of brevity, we will not disuss this 
here. 



IV Lattice Input for CKM Fits 

Table ^ shows the input from the lattice, experiment and elsewhere used by 
us and compare it with the works of Ciuchini et al |S2| and Hocker et al 
|33) . The corresponding determination of the CKM parameter p, fj and unitarity 
angles a, /?, 7 as well as several other quantities of interest are also shown. Note 
that our error on fBdVSBd and on ^ are appreciably bigger than used in the 
other studies. This is especially so for ^, where for quite sometime we have been 
cautioning that the error of w 0.05 that was commonly taken was a serious 
underestimate Recently Kronfeld and Ryan and Yamada jJB] have 

also argued for a reassessment of errors on ^ due to the presence of chiral logs. 
Following this development as of LAT'02 larger error on ^ is now being widely 
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Table 1: Comparison of some fits. 



Input Quantity 


Atwood & Soni 31 


Ciuchini et al 32 


Hocker et aV'i'i 


/in . — - , / l/^ Ji 
'^UC 1 * 110 1 ' CO 1 


.085 ± .017 


.089 ± .009 


.087 ± .006 ± .014 


Fb,\Ibb, MeV 


230 ± 50 


230 ± 25 ± 20 


230 ± 28 ± 28 




1.16 ± .08 


1.14 ± .04 ± .05 


1.16 ± .03 ± .05 


Bk 


.86 ±0.15 


.87 ±0.06 ±0.13 


.87± .06± .13 


Output Quantity 








sin 2/3 


.70± .10 


.695 ± .065 


.68± .18 


sin 2a 


-.50± .32 


-.425 ± .220 




7 


46.2° ±9.1° 


54.85 ±6.0 


56± 19 


»7 


.30± .05 


.316 ± .040 


.34± .12 


P 


.25± .07 


.22 ± .038 


.22± .14 


\VtdlVts\ 


.185 ± .015 




.19± .04 




19.8 ±3.5 


^ ' ■•J-0.7 


24.6 ±9.1 


Jcp 

BR{K+ -K+uD) 
BR{Kl TT°VV) 


(2.55 ± .35) X 10-5 
(0.67 ±0.10) X 10-10 
(0.225 ±0.065) X lO-^" 




(2.8 ± .8) X 10-5 
(.74± .23) X 10-1" 
(.27± .14) X 10-1" 



advocated; for example, Lellouch ^| in his review at ICHEP02 summarized 
^ = 1.18±0.04ti^ 

The SM fits now give (sm2/3)'5A^ = 0.70 ± 0.10 as weU as allowed ranges for 
7,77, p etc (see Table While these fits provide fairly restrictive range for (3 
and 7, a is constrained rather poorly. Note also that Bg — B^ mass difference 
AtoBs = 19.8 ± 3.5ps"i is now constrained with a one sigma accuracy of about 
15%; measurements at the Tevatron and later at the LHC should be able to 
test this important prediction of the SM. Meanwhile measurements of the CP 
asymmetry in i? — > tpK^ is already providing quite an impressive determination, 
sin 2/? — 0.734 ± 0.054 [37] in good agreement with the theoretical prediction. 
It is important to note also that just in the past year B-factory experiments 
have improved the determination of sin 2/3 from an error of ± 0.10 down to 
± 0.05. With the anticipated increase in luminosities of the B-factories, along 
with results from the Tevatron, this error should go down further to « 0.02 in 
another year or two. (Recall that the intrinsic theory error in the determination 
of sin 2/3 is expected to be about < 0.01) 38 

It is instructive to reflect on the pace of theoretical progress in constraining 
sin 2/3'^*^. For this purpose we may compare the inputs used in fits of ~ 1995 
|39|. with that of « 2001 PJ. Indeed the 2001 fit has reduced the error on 
sin 2/? from 0.20 to 0.10; correspondingly the error on fj and on Ato^^ is also 
appreciably reduced. However, only some of this improvement can be related 
directly to lattice computations. In fact it seems that a large portion of the 
improvement is due actually to the reduction especially in the error on Vcb and 
to some degree on Vub/Vcb wherein the role of the lattice is less clear. 
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What should be clear is that it will be extremely difficult to reduce the 
theory error on (sin 2(i)^^ from the current level of ± 0.10 down to level of 0.02 
that the experiment is anticipated to reach in the very near future. Thus to 
test the SM more precisely will require clean determination of the other angles 
a and 7 directly from experiment. We will come back to this point later on. 

IV. 1 Important Input from Experiment on the horizon for 
CKM Determination 

1. S± ^ T± + Vr{i>r) 

With 10* or more B - B pairs that BELLE and BABAR each will soon 
have access to, an experimental determination of (actually Jb x V^b) 
may be feasible. Using from the lattice Jb ~ 207 ± 35Af eV^^ and 
Vub/Vcb « 0.085 ± 0.017 one gets an estimate, Br [B ^ t + Vr ~ (7.8 ± 
2.0) X 10~^). Decays of r into final states with {v' s)+ ji (c, p or tt) have 
a total branching ratio of around 50%. So with a few percent detection 
efficiency there should be a few hundreds of events for B^ + Vrivr), 

a respectable sample to provide a reasonable determination of Jb x Vub 
and an important check on the lattice calculation. 

2. B Iv-y 

Unlike B —^ Iv, Ivj {I — e, /i) does not suffer from helicity suppression 
although it is suppressed by a. Emission of the photon from the light 
quark also tends to enhance the process although precise calculation of 
the Br is difficult to make 001) estimates are in the range of 1 — 6 x 
10~^, i.e. about an order of magnitude more than the 2-body helicity 
suppressed modes, B —f Iv. The constituent quark model, although too 
simple to provide reliable details, perhaps does give a valid qualitative 
picture indicating a "hard" photon spectrum |40| : 

where = Ej/tub, and yields a total Br w 5 x 10~^ with a constituent 
light quark mass of about 350 MeV and /b = 200Afey. Predictions 
from several other estimates are given in Table |2 These radiative modes 
should be accessible accompanied by /i or e with the data samples currently 
available. In making contact with the phenomenological models, energy 
spectra of the photon and of the neutrino [ i. e. the invarriant mass of (7+ 
the charged lepton)] would be especially useful. Detailed experimental 
studies of these radiative decays would also give another handle on the 
approximate value of /b- 



3. B^ p + -f 
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Table 2: Sample of rough estimates for Ivj. 



Method 


Reference 


Br 


const, quark model 
hght cone + HQET 
light-front 
Rel. potential 
QCD-fact,SCET 


Atwood et al 40) 
Korchemsky et al|41) 
Ceng et algT] 
Colangelo et al'H' 
Descotes-Genon and SachrajdaEH 


w 5 X 10"" 
w 5 X 10-6 
w 4 X 10-6 
w 1 X 10-6 
« afew X 10^6 



Another important input from experiment that could aid in the determina- 
tion of the CKM-parameters (esp. | Vtd\/\ Vts\) is the reaction ^ p + 7 
since 

Br{B° ^ p°7) ^ [1 - ml/ml] 
Br{B K*j) [1 - m|.. /m|] 

The expected Br{B^ P + l) ~ 1 x 10^^ seems within reach of ex- 
periment. An accurate calculation of the (SU(3)) breaking ratio of form 
factors [rf-*''(0)]/[Tf-"-^'*(0)], for example by lattice methods, the fea- 
sibility of which was demonstrated long time already 02 QHI > along with 
the anticipated experimental measurement of _B — > p -I- 7 could lead to 
another determination of Vtrf /Vt» 

This method of extracting Vtrf/Vts has some advantages and some dis- 
advantages compared to Bg — Bg oscillations. For B ^ the relevant 
operator is a bilinear one, whereas for Bg — B^ oscillations a 4-quark oper- 
ator enters; the renormalization of a lattice 4-quark operator can be more 
complicated compared to a bilinear one. On the other hand, _B — )■ p + 7 
involves a large recoil thereby extracting the form factor at or near ~ 
(where q is the 4-momentum of the photon) is numerically difficult on the 
lattice. 

From the experimental side in the numerator of eq. (|13|l only neutral 
B's, B'^,B^ — > p'^7 should be used as charged B^ p^-f provides a 
non-negligible long-distance contribution 03] which is proportional to Vub 
(and independent of Vtd) coming from the annihilation graph and cannot 
be estimated accurately. 

Due to these anticipated input from experiments along with developments 
in theory, especially with the expected improvement in computational resources 
because of the Scientific Discovery Through Advanced Computing (SCIDAC) 
initiative (5^, it is fairly safe to say that errors on SM parameters such as sin 2/3 
will go down by a factor of about 2-3. However, with larger pool of B-samples 
that are expected from B-factories and the hadron facilities, experiments should 
be able to directly determine sin 2/3 from CP asymmetry measurements in B — > 
ijjKs to an accuracy of 0.02; so experiment is likely to stay ahead of theory. 



(0) 



Vtd 



(13) 
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Uncovering new sources of CP violation in B-physics though may well require 
more precise tests, as we will emphasize in the next few paras. 

V Theory Errors and the hunt for new sources 
of CP violation in B-physics 

Theory errors should be a concern as they can thwart experimental efforts to 
search for the beyond the standard model (BSM)-CP odd new phase(s) which 
we will collectively denote as %. The main point is that x "^^Y 'well cause only 
small deviations from the SM in B-physics. Indeed the emerging understanding 
of the CKM-paradigm serves as an important lesson in this regard. Thc-CP-odd 
phase (5i3, in the standard notation, ((5i3 « 7 « 50 ± lOdeg) is not small and 
although it causes 0(1) CP-asymmetry in i? ^ 'ipKs its effects in K]^ tttt, e or 
e' are very small, O{10~^ — 10~^) respectively [33]. Analogously it is clearly not 
inconceivable that even though a BSM-CP-odd phase x is not small its effect on 
B-physics will be small. As an example, this may happen if x arises in models 
with extra Higgs bosons; then its effects may be much larger in top physics and 
quite small in B-physics [T7j . 

For one thing this means we may need very large data samples of B's. Indeed 
for an asymmetry of O(10~^) (as in Kl), since the relevant Br is unlikely to be 
larger than « 10"'^, which is about the branching ratio for B v'-^s, detection 
may require O(10^°) B's. Higher luminosity super BELLE/BABAR B-factories 
as well as efforts at hadron B-facilities BTEV and LHCB may well be needed 
in the hunt for x- 

In the search for x the ability to detect small deviations from the SM also 
though requires that we develop tests of the SM that use little or no theory 
assumptions and are as free of theory errors as possible. Note in this regard 
that for detection of deviations from the SM at the level of « 10~^ means that 
even isospin approximation, widely used in many methods for extracting angles 
of the unitarity triangle can mask x a-nd thereby defeat the experimental effort 
for detection of new physics. 

Motivated by these considerations wc now discuss methods for getting uni- 
tarity angles with very little theory error, i.e. to 0{< 1%). 

V.l Pristine Determination of the Unitarity Triangle via 

B -> KD^ 

Angles of the unitarity triangle can be obtained very "cleanly" i.e. without 
any theoretical assumptions from analysis of final states containing D'^,D^ in 
charged or neutral B-decays|2|. 

7 can be extracted from a study of direct CP violation in charged B-decays, 
B± ^ K^D°, Z)o mi EH EOI. (5 = (/3 - a -t- tt) = 2/3 -h 7 ISH 112 E3 EH as 
well as (3 can be obtained from time dependent CP-asymmetry measurements 
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in B°,B° K^D°,D°. In both cases common final states of D°,D° have to 
be used, as flavor tagging of D°,D^ is very difficult^^l- There are 3 types of 
such common final states: 

1. D'^,D'^ decays to CP-non-eigenstatcs that are doubly Cabibbo suppressed 
|49| . for example, K~^Tr~, K'^p^, K*^tt~ etc. 

2. D^,D° decays to CP-eigenstates gHl, for example, K+R-, tt+tt", KItt^ 

3. D'^,D'^ decays to CP-non-eigenstatcs that are singly Cabibbo suppressed 
120] for example, K*+K^, K^K*, p'^n^ , p^ir^ etc. 

It turns out that CP asymmetry are expected to be small (< 10%) for CPES 
and for singly Cabibbo suppressed modes (i.e. 2nd and 3rd type) whereas the 
interference and CP asymmetry is maximal for CPNES (1st type). On the other 
hand, the branching ratios are expected to be largest for CPES and smallest 
for doubly Cabibbo suppressed CPNES modes. The general expectations are 
that for extraction of 7, doubly Cabibbo suppressed modes should be most 
eflicient among the three types. However, this is not guaranteed and all three 
methods should, for sure, be used. What is important, for the long run is to 
note that only two common decay modes of Z3°Z)° are needed to give enough 
observables to algebraically solve for the CP-odd weak phase 7, the strong 
final states phase(s) as well as the suppressed Br(i?~ K~D^) that is very 
difficult to measure experimentally; indeed perhaps a dozen or so such modes 
are available. This should greatly help the analysis in extracting a precise value 
for 7 without discrete ambiguities. 

For time dendent CP asymmetry ^4 in ,B'^ K'-*D'-\ Z)°, the discussion 
is analogous to the above. Again for extraction of S (as well as (3) one needs only 
two common final states of 13°, from the many; whether they be CPNES, 
CPES doubly or singly Cabibbo suppressed modes, that are available; however, 
both modes cannot be CPES. 

Especially noteworthy is the fact that for clean extraction of the angles, 
final states containing D^,D^ in the decays of B^,B^^B^ are involved and 
furthermore common final states of , 13° decay play a critical role and should 
aid in increasing the efficiency of the analysis. 

In passing we briefiy note of the analogous methods involving Bs decays 
to Ds (or their vector counter-parts|56j') that can also give 7 very 

cleanly [57]. 

VI Exact Chiral Symmetry on the Lattice 
VI. 1 Introduction 

In the past few years a significant development for lattice gauge computations 
has taken place. For the first time, we have practically viable discretization 
methods that exhibit exact chiral symmetry on the lattice even at a finite lat- 
tice spacing, i.e. even before the continuum limit is taken. By now, not only 
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the viability of these methods has been convincingly demonstrated, large scale 
simulations, with some success, have already been using them to address some 
outstanding problem in weak interaction phenomenology pertaining to K 2tt 
that were very difficult to address heretofore, as briefly reminded below. 

VI. 2 Difficulties of calculating weak matrix elements on 
the lattice with conventional discretizations 

Recall that conventionally there are two fermion discretizations: Wilson and 
Staggered (or Kogut-Susskind). Wilson fermion explicitly break chiral symme- 
try whereas conventional staggered fermions while possessing some residual chi- 
ral symmetry break flavor symmetry. The difficulties of maintaining chiral sym- 
metry on the lattice is enunciated in the form of a no-go fNielsen-Ninomiva|58| ) 
theorem. 

Although conventional wisdom says that these symmetries get restored on 
the lattice in the continuum limit, in practice, in the study of hadronic weak 
decays, lack of chiral symmetry imposes an extremely serious if not an insur- 
mountable limitation. 

Lack of chiral symmetry leads to two types of significant difficulties: 

1) Precise renormalization of 4-quark operations can become a difficult fine- 
tuning problem. The point is that, in the absence of chiral symmetry operators 
such as Oll = [s7At(l — 75)c^]^, that are relevant to K — K oscillations and Bk 
computation, mix under renormalization with wrong chirality operators 159] , for 
example, with Opp = (575^)^. 

The problem is that whereas < K\Oll\K > is proportional to the quark 
mass and therefore vanishes in the chiral limit, < K\Opp\K > goes to a constant 
in the chiral limit. Thus even if the mixing coefficients of the wrong chirally 
operators are small you need to know them very accurately in order to precisely 
extract the matrix element of physical interest. 

2) Mixing with lower dimensional operators is even a worse problem. This 
happens, for instance, when one considers the operators of the AS* = 1 Hamil- 
tonian (e.g. S7^(l — 75)um7^(1 — 75)^) relevant to K ^ 27r. Now such a dim- 
6 operation mixes with lower dimensional operations, for example, sd, sj^d, 
,S(7y,„(i(7^'^|(j[)|. The mixing coefficients are now power divergent, for example 
~ a~" (n=3 for sd and 575^ and n = 1 for sa^i^dC^"). So they become in- 
creasingly important in the continuum limit. Non-perturbative methods (that 
respect chiral symmetry) are essential for handling them. This was the main 
reason that early efforts (BT] to calculate A/ = ^, K ^ 2tt amplitudes on the 
lattice did not make much progress. 

VI. 3 Domain Wall Fermions 

In 1992, Kaplan j221 in a celebrated paper showed a simple method to attain 
exact chiral symmetry on the lattice even at finite lattice spacing. This remark- 
able feat is accomplished by embedding the 4-dim theory on to 5-dim with a 
fermion mass-term that has the shape of a domain wall across the 4-dimensional 
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boundary and switches sign. The low lying mode bound to the walls then pos- 
sesses exact chiral symmetry on the lattice in the limit that the length (ig) 
of the 5th dimension has an infinite extent Nielsen-Ninomiya theorem 
is evaded as in Kaplan's construction the number of fermionic degrees of free- 
dom per (4-dim) site is no longer finite as assumed in the theorem, and in fact 
becomes infinite as Ls — > oo. Narayanan and Neuberger jM] gave an elegant 
flavor interpretation of this fifth dimension. 

In 1997 1st (quenched) lattice QCD simulation done to test the practical 
viability of this approach showed very encouraging results . In those numer- 
ical simulations for QCD actually the domain wall formulation of Shamir |ti3 | 
was used. These early results showed that even with a modest extent of the 5th 
dimension, domain wall fermion possesses very good chiral behavior setting the 
stage for their use in large scale simulations. 

Since DWF are continuum like their renormalization (perturbative and non- 
perturbative) properties are fairly simple. Also discretization errors tend to go 
as O(a^) rendering them with very good scaling properties which tends to offset 
the cost of the extra-dimension. 

Since in practice the extent of the 5th dimension is finite, the coupling be- 
tween the two walls separated by Lg causes a coupling between the light modes 
and gives them a residual mass, rrires- This mass can be measured quite pre- 
cisely inSl ESI . In low energy applications one can systematically include the 
effect of rrires in the context of an effective chiral lagrangian |b7| . 

VI. 4 Application to Bk 

CP-PACS 22 and RBC221 collaborations have made considerable progress 
towards a precise calculation of Bk with DWF. Both results are in the range 
of l-2a below the old result from JLQCD "T, Bk[2 GeV]=0.628 ± 0.042. CP- 
PACS and RBC central values for Bk differ by about 5-10%; most likely this 
difference is due to the fact that CP-PACS uses 1-loop lattice perturbation 
168" "HS* theory for renormalization of the AS* = 2 operators whereas RBC is 
using non-perturbative renormalization |67| . Efforts are now underway to repeat 
this calculation at weaker (quenched) coupling [7D| as well as with dynamical 
domain wall quarks |15) . 

VII AI =1/2 Rule and e'/e: Progress and Out- 
look. 

VII. 1 Introduction 

There have been two recent attempts by the CP-PACS |2S1 and the RBC [E] 
collaborations, at attacking this old problem on the lattice using the relatively 
new discretization method of domain wall fermions (DWF) |62l 1631 

First lattice studies of K tttt by both CP-PACS |2S| and RBC [IS| with 

DWF used the lowest order chiral perturbation theory (ChPT) approach sug- 
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gested by Bernard et al. [7J|. The method then calls for using the lattice to 
compute matrix element of 4-quark operators between K tt and K —^ vacuum 
which are used to obtain the corresponding desired — > 27r matrix elements |72). 
While this is a simple method which avoids technical (Maiani- Testa Theorem 
[73,) and also practical, computational limitations, it is nevertheless a severe 
approximation. In particular, at the leading order in ChPT being used, final 
state interactions, which in reality are very likely important j74| are necessarily 
absent. Note also that in this approximation the good chiral behavior of domain 
wall fermions becomes crucial. For one thing in the absence of chiral symme- 
try, the unphysical (cubically divergent) contribution from mixing with lower 
dimensional operators to {tt\9^-'^\K) cannot be subtracted away in a relatively 
simple way by using {0\9^'-^\K°). Furthermore, the renormalization of 4-quark 
operators also becomes vastly more complicated due to the mixing of operators 
with the wrong chirality ones|59|. Since the K —> tt matrix elements of some 
of the operators go to a constant in the chiral limit, whereas those of the right 
chirality tend to vanish, subtraction of the unwanted contribution needs to be 
done at a very high precision, i.e. it becomes a fine tuning problem. For these 
reasons, as mentioned above, earlier efforts jM] for computation of if — > 27r and 
e'/e to the LO in ChPT by the use oi K ^ tt and iiT — > on the lattice with 
Wilson fermions, which explicitly break chiral symmetry, were able to make 
little headway. 

Since chiral symmetry is so critical in the calculation of matrix elements for 
K 277 and since, due to the finite extent of the 5th dimension, rigorously 
speaking, domain wall quarks do not possess exact chiral symmetry, it is im- 
portant to be able to take into account residual symmetry breaking effects in a 
systematic fashion. For matrix element dominated by long-distance physics this 
can be accomplished by shifting the bare masses in ChPT by rrires , where m^es 
is the residual quark mass which the massless quarks on the lattice possess due 
to the coupling between the walls of the 5th dimension [ST]- On the lattice we 
can calculate rrires quite precisely |65[ I66| and the chiral limit is then taken by 
setting {niguark + rrires) ~^ 0. 

For operators such as which receive power divergent (i.e. short distance) 
contributions that are not physical and have to be subtracted away, the sym- 
metry breaking effect cannot be precisely described by rrires |13| . Fortunately, 
the LEG can still be computed accurately by taking the slope of the matrix 
elements with the rriquark so long as rrires is independent of rriquark to a good 
approximation jl3| : of course this does require that in actual simulations the 
length of the 5th dimension is sufhciently long that rrires << rriquark- 

Note also that for power divergent subtractions ChPT is taken into account 
to all orders I^ . 

Table 13 gives the (subtracted) K ^ tt matrix clement of all of the 4-quark 
operators of interest for the I — 1/2 and the 7 = 3/2 channel. This can be used 
to obtain the K ^ 2Tr amplitudes via eq. (201) of Ref. |13| . 

Table H gives the full resuhs for ReAo, Re^2, w"^ = and e'/e of 

RBcna. 

While both the groups [2211131 use LOChPT, DWF and the quenched approx- 
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Table 3: The lattice values for the low energy, chiral perturbation theory con- 
stants decomposed by isospin for Qi to Qiq. (Taken from |13| ) 



i 


"i,lat 


"i,lat 


1 


-1.19(31) X 10"^ 


-1.38(6) X 10-« 


2 


2.22(16) X lO-'^ 


-1.38(6) X 10-6 


3 


0.15(113) X 10-5 


0.0 


4 


3.55(96) X 10-5 


0.0 


5 


-2.97(100) X 10-5 


0.0 


6 


-8.12(98) X 10-5 


0.0 


7 


-3.22(16) X 10-6 


-1.61(8) X 10-6 


8 


-9.92(54) X 10-6 


-4.96(27) X 10-6 


9 


-1.85(16) X 10-5 


-2.07(9) X 10-6 


10 


1.55(31) X 10-5 


-2.07(9) X 10-6 



Table 4: Final values for physical quantities using 1-loop full QCD extrapola- 
tions to the physical kaon mass and a value of fj, — 2.13 GeV for the matching 
between the lattice and continuum; Taken from 13 . (Statistical errors only ) 







This calculation 


Quantity 


Experiment 


(statistical errors only) 


Reylo(GeV) 


3.33 X 10- '' 


(2.96 ±0.17) X 10- ^ 


ReA2(GeV 


1.50 X 10-s 


(1.172 ± 0.053) X 10-8 




22.2 


(25.3 ± 1.8) 


Rc(e'/e) 


(15.3 ±2.6) X 10--^ (NA48) 


(-4.0 ±2.3) X lO-'' 




(20.7 ±2.8) X IQ-'' (KTEV) 





imation, there are some important differences in these two calculations as well. 
For one thing, RBC 13 used the standard Wilson gauge actions whereas CP- 
PACS I^S, used renormalization group improved (Iwasaki) gauge action |75| . 
Also in their extractions of Re^2 (and Bk), RBC used the 1-loop quenched 
chiral perturbation theory ^7S^ to fit (TrlQ^^^lJ-iT) whereas CP-PACS used a phe- 
nomenological fit. 

VII.2 When the Dust Settles 

I. Regarding Re e'/e 

a) Key Contributions. 

Listed below are the key contributions to e'/e from / = and 2 
final states, all given in miits of 10—* resulting from^l^j. Recall, 
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experiment finds (in this unit) [771 [7H| Re e'/e = 17 ± 2. (Note 
contributions not shown are neghgible in comparison). 



Operator 



/ = 



/ = 2 



Qa 



-4.8± 1.1 
14.2 ± 1.9 
1.48 ± .12 



-16.97± .84 
1.56 ± .00 



Qe and Qs are not the only ones that matter. 
Although, as widely expected (75]-|S2l7 Qe and Qs are the dominant 
players, due to the cancellations between these two contributions, 
other operators (e.g. Q4) seem to be making an appreciable difference 
to the final result. 

Buras approximate formula |81| . 
In this context recall Buras' approximate formula 



From our lattice data one can see that the contribution of operators 
other than Qg and Qs is about 60% of e'/eje+s- 
Cancellation not between large numbers. 
While there is a cancellation between Qq and Qs , in magnitude each 
of this contribution is comparable to the experimental number for 
e'/e. Had it been that 



then the cancellation would have been between "large numbers" (com- 
pared to the final result that one is seeking) and the prognosis for 
future improvements would have been even harder. 

Unnatural Cancellations. 
The substantial cancellation (~ 85%) between contributions of Qq 
and Qs to e'/e, in all likelihood, is not natural, i.e. not stable to per- 
turbations. Recall that these numbers emerge after using at least 3 
key (uncontrolled) approximations: lowest order chiral perturbation 
theory, quench approximation and heavy charm quark. It is virtually 
impossible that these approximations affect Qe and Qs in the same 
way. Indeed there are good reasons to think that both chiral pertur- 
bation theory and quench approximation are having a bigger effect 
on Qq than on Qs- It seems reasonable therefore to expect that in 
improved calculations of e'/ e these cancellations between Qq and Qs 
will not remain. 



e'/e^e'/e\Q+s^e'/e\Q + e'/e\s 



(14) 



e'/ele, e /e\s » e /e\expt- 



(15) 
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f) Phenomenological bound. 

To the extent that e'/£|/=2 < 0, there is a useful pheomenological 
bound, 

e7e|/=o > e'/elertpt (16) 

with which one can test the SM. For the purpose of this test the 
cancellations between the / = and 1 = 2 contributions are not 
quite relevant. Note that our current data gives left hand side of 
eqn.(16) of '--^ 11 ± 2 whereas the RHS (from experiment) is ~ 17± 2. 
Clearly then e'/e|/=o rnust increase appreciably as improvements in 
our lattice calculation are made if the SM's description of CP is to 
continue to hold. 

These considerations suggest that tests of the SM with improvements 
in accuracy appear feasible. 

Repercussions for the origin of the A/ =1/2 Rule. 

The octet enhancement, i.e. §^4^ ^ 20 >> 1, has been a long stand- 
ing puzzle in Particle Physics. The lattice calculation with domain wall 
quarks, although not having sufficient control over all the systematic er- 
rors cannot at present give a reliable result for e'/e, they do provide with a 
useful and unambiguous information for the A/ =1/2 enhancement. The 
lattice result leads to an important and remarkable conclusion regarding 
the A/ =1/2 rule as can be seen from Table The contribution of 
ReA2 and especially of ReAo originate almost entirely from the aboriginal 
4-fermi operator Q2. Indeed, we find 

Operator Re^o (GeV) ReAa (GeV) 

Qi (3.48 ± .77) X 10-8 (-.363 ± .016) x IQ-^ 

Qa (24.5 ± 1.6) x 10"^ (1.520 ± .068) x 10"** 

Qe (0.050 ± 0.006) X 10-8 

These numbers should be compared to the experimental ones: HeA^^^^ = 
33.30 X 10-8 GeV, Re^2"'^* = 1.50 x 10"® GeV. Clearly ReAo is completely 
dominated by Q2, making about ~ 80-85% contribution and Qi makes 
the remaining ~ 15% contribution to HbAq. In particular the contribution 
of Qe to ReAo is completely negligible, being ~ 0.2%. This is in sharp 
contrast to some model calculations for e'/e in which contribution of Qe to 
Re^o is typically almost 20-30% jZHI IHOI ■ Note that while numbers given 
here for contribution of individual operators are based on calculations at 
/i ~ 2 GeV we have studied the /i dependence from ~ 1.3-2 GeV and the 
dependence is quite mild |T5] . 

It must be emphasized that all these quantitative findings, in particular Qe 
vs Q2 contributions to K ^ 2tt{I = 0), are based on LOChPT calculations 
and this could change as higher order corrections arc included. 
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It is indeed interesting and ironical that the penguin operators |S3| origi- 
naUy invoked to explain the A J =1/2 rule seem to play little role therein, 
at least in the context of our lattice calculation^Jl- However, the subse- 
quent conjecture of their importance to rendering a largish e'/e seems 
to be substantiated although significant theoretical progress still needs to 
be made before the repercussions of the precise experimental measurement 
can be fully assessed. 

VII. 3 Approximations and Concerns 

There were several approximations made in the lattice of calculations [231 [T5| 
using DWF that were recently completed. 

1. The quenched approximation so that quark antiquark loops are ignored 
in the propagation of the gauge field (gluon). 

2. Lowest order chiral perturbation theory (LOChPT), so that the matrix 
elements of the 4-quark operators are calculated in the leading order (LO) 
in this approximation. This means that operators [Qi-Qe, Q9, Qw] that 
transforms as (8,1) and/or (27,1) under SU(3)LxSU(3)i^ are calculated to 
0{p'^) whereas those of Q7, Qs which transform as (8,8) are to 0(p° — > 1) 
in the chiral expansion. 

3. The charm quark is assumed to be very heavy and integrated out. The ef- 
fective Hamiltonian [HHIIHEIj consequently consists of only 3 active flavors: 

M, d, s. 

These approximations are uncontrolled, i.e. we do not have a reliable es- 
timate of how inaccurate they are. While the quenched approximation seems 
to be accurate to 10-15% in many spectrums and decay constant calculations, 
it may well be a lot worse for some hadronic matrix elements. In particular, 
comparison of the analytical formulas for {TT\Q^^^g\K) to NLO ^87 in full ChPT 
with the lattice data 13 obtained using the quenched approximation shows 
that the logs of the full ChPT seem to be absent Also the matrix element 
of Qq, which is of crucial importance to e'/e is claimed to be very susceptible 
to quenching effects 

Although, in many low energy applications, LOChPT works fairly well, in 
K ^ 2Tr there are reasons to be suspicious. First of all an important mass 
scale here is ttik and not just m-^. Furthermore, in the 1 = channel the tt-tt 
rescattering effects (FSI), which cannot occur in the LOChPT, are likely to be 
quite important [Tl]- Indeed for {ttttIQqIK) higher order chiral corrections may 
well be intertwined with a O^'^(cr) resonance in the tt-tt channel [101 ED- 

Since the mass of the charm quark is only ~ 1.3 GeV, integrating it out as- 
suming it is very heavy (i.e. >> Aqcd) is very likely not a good approximation. 
Corrections from higher dimensional operators are likely to be sizeable ^2^. Also 
in the 4-flavor theory, GIM cancellation forbids power-divergent mixing of dim-6 
operators of the / = 0, -ffe// with lower dim operators, so the 4 flavor theory is 
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preferable over the 3 flavor theory for that reason too, although this advantage 
is only relevant in computation of KeAo In the calculation of Im^o 

where (8,1) penguin operators, such as Qq, become important which originate 
from integrating out the top quark, the top quark is so heavy compared to the 
lattice cut-off that integrating it out is unavoidable. 

VII.4 Future Outlook 

While the above approximations used in the current lattice calculations are not 
controllable, systematic improvements are feasible and efforts at these are well 
underway. 

First, recent studies of renormalization group-improved gauge actions in the 
context of domain wall fermions appear very promising in significant improve- 
ments in chiral symmetry (which was already remarkably good) ||96j. Efforts 
are also underway towards creating a large ensemble of gauge configurations 
with dynamical (2 flavor) domain wall quarks jl6| with lattice spacing ^ (2 
GeV"""^). This should allow a first study of the quenching effects in X 27r 
matrix elements in another year or so. Also work is being done at finer lattice 
spacing ~ (3 GeV~^) with the hope that this will allow a better treatment of 
the charm-quark and a calculation of the matrix elements in the effective theory 
with 4 active flavors {u, d, s, c) |70j . 

Note also that new calculations of the K ^ 2tt matrix elements have be- 
gun 1^ using another discretization (overlap fermions ,98j) possessing excellent 
chiral symmetry. 

Recent works also show how lattice computations of all the matrix elements 
relevant to X — s- 27r and e'/e can be obtained beyond the leading order in 
ChPT. In one method 88 matrix elements of all the relevant operators (A/ = 
1/2 or 3/2) can be obtained to NLO by using lattice computations of K-K, 
K TT, K —f and K —> 2tt ai the two unphysical kinematics (tok — 
and rriK = 2711^^) wherein Maiani- Testa Theorem 73 can be evaded. In another 
construction |99j the K —> 27: matrix elements for A/ — | transitions can be 
obtained to NLO by using lattice computations of K — )■ 27r with momentum 
insertion on one of the final state pions. 

Indeed in a very interesting paper Lellouch and Luscher have also proposed 
a method wherein K ^ 2tt matrix elements may be directly calculated with- 
out using ChPT by relating them to finite volume correlation functions jlUU) . 
There is also a proposed method which makes use of dispersion relations to cal- 
culate physical K tttt amplitudes to all orders in ChPT [T7ITj . We note that 
both of these methods make rather stringent demands on unitarity therefore 
their implementation, especially for the A/ =1/2 case, may need full QCD 
simulations. 
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VIII Unitarity Triangle from A'-Decays 



While determination of the unitarity triangle (UT) from B decays has been 
receiving considerable attention and is much in the news, it is useful to note 
that not only an independent determination of the UT can also be made purely 
from if-decays but it is important to do so. In principle, there are four physical 
processes that can be useful here, whereas any three of them would be sufficient. 

1. e, the indirect CP- violation parameter characterizing the CP- violation in 
Kl 27r. 

2. The Br[K^ — )■ tt+j/S'), a determination of which has been underway for a 
long time and a crude measurement now exists thanks to the two candidate 
events that have been seen, Br{K+ n+vv) = (1.57t^g2^) x 10"^° [102. 

3. There is considerable experimental interest in measuring the Br^K^ 
■K^vv). This is a very interesting mode which is CP violating |103j and 
is theoretically extremely clean, but clearly an extremely difficult experi- 
mental challenge. 

4. The direct CP-violation parameter e'/e. Although the experimental num- 
ber is now quite precisely known [771 178| . it can only be useful in the 
context of the UT, if the theory can be brought under-control. Renewed 
interest on the lattice, in light of recent progress in maintenance of chiral 
symmetry on the lattice (described briefly in preceding pages) gives one 
some encouragement that perhaps a few years down the road we would 
be able to make use of the experimental result and translate it into the 
CP violation parameter r] of the CKM paradigm. In the p-rj plane, e'/e, 
when the numerical value of e is taken from experiment, would provide a 
horizontal line (actually a band due to the error in theory and in experi- 
ment). 

As is well known transplanting e to p, 77 plane does require knowledge 
of the non-perturbative hadronic parameter, Bk- Fortunately as already 
mentioned, lattice calculations oi Bk are now quite mature. In fact several 
different discretization methods have been used to determine this impor- 
tant quantity. While the current accuracy is around 15%, efforts with 
dynamical quarks are underway and in 3-5 years we should expect the 
accuracy to improve appreciably. 

The theory for K'^ tt^i'D is also rather clean The basic process 

s — *■ dviy is dominated by the top quark. Conversion of that to —^ ■n'^vv can 
be done using isospin by relating it to the thoroughly studied charge current 
process K — > -nev. 

The observed Br, deduced on the basis of the 2 events seen so far, is con- 
sistent with the expectation from the SM; our |2] CKM fits give, BR{K^ 
TT^iyp) = (0.67 ± 0.10) X 10~^". Experimental efforts are underway to improve 
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this measurement in the near future at BNL and further down the road at FNAL 

The decay n^vv is fascinating as it is CP violating. The the- 

ory 121] in this ease is even clearer then for the charged counter part and 
BR{Kl K°vv) = 1.5 X lO-^A^AiOTy^ Qur |3T] fit value is BR{Kl 
TT^vv) = (.23 ± .07) X 10^^°. So an experimental measurement would give a 
clean determination of the CKM phase 77. Note though that {A « Vcb) the 
current accuracy in A(~ 7%) should be improved otherwise it introduces signif- 
icant error on the r] determination. The Kl n^vv experiment is clearly very 
challenging and it is receiving attention at KEK (E391), at BNL (KOPIO) and 
at FNAL (CKM)il03i- 

Given the intrinsic difficulties of this experiment and those of an accurate 
theoretical calculation of e'/e it would be interesting to see which of these is 
brought under control first. 
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